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The molecular structure of a 3d 5 configuration ion in a trigonal ligand field is theoretically established on the basis 
of the 252x252 complete energy matrice. The optical absorption and the electron paramagnetic resonance spectra of 
the Fe i+ ion in the sapphire crystal (a— AI2O3) have been studied by diagonalizing the complete energy matrice. The 
calculated results are in very good agreement with previous experimental observations. The strength of the transition 
probabilities between pairs of the energy levels have been calculated to determine the possibility to achieve a population 
inversion in the ground state by applying optical pumping to the crystal. Preliminary results based on the computed 
transition parameters and on the maser rate equations model show that a 31 GHz maser signal can be effectively 
generated depending on the cryogenic resonator design. 
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I. INTRODUCTION 

Masers based on paramagnetic ions hosted in a crystal 
were extensively studied and used as stable microwave 
sources or as low noise amplifiers in the 60's . Maser action 
was obtained at low temperature in different doped crystals, 
which the most popular was the sapphire doped with Cr 3+ 
or Rubis. Numerous types of maser were designed generally 
based on a doped crystal inserted into a waveguide or a 
metallic cavity cooled down the liquid Helium temperature. 
The doping level was typically of the order of 0.05% or so 
and for a X-band source the cryogenic resonator presented a 
Q-factor of the order of 10,000 or less. Most of the systems 
were based on the ground state levels multiplicity to get 
the required population inversion through the action of a 
microwave pumping. Some attempts to optically pump the 
crystal have been realized^. But no real advantage compared 
to the classical maser scheme has been demonstrated. The 
lack of narrowband and easy to handle laser source as well 
as the poor quality of available doped crystals were certainly 
the main difficulties to achieve a reliable optically pumped 
microwave maser. Eventually the interest for this technology 
was almost stopped after the development of low noise figure 
cryogenic semiconductor amplifiers, stable quartz crystal 
oscillators and low noise frequency synthesis. 

More recently, we reported the possibility to get a 12 GHz 
maser signal presenting a short term relative frequency in- 
stability below 1 x 10~ 14 at short term with a large marge 
of progress^. Such a frequency stability performance is re- 
quired for a number of very demanding applications as metrol- 
ogy, space navigation, radioastronomy or fundamental physic 
experiments. The search for a reliable microwave source 
presenting a relative frequency stability of some 10~ 15 is a 
challenge that motivates a lot of innovative works all around 
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the world*—. In our first prototype the sustaining amplifica- 
tion is achieved in a cryogenic sapphire whispering gallery 
mode resonator containing a small amount of paramagnetic 
Fe 3+ ions. The Fe 3+ ion exhibits three ground state energy 
levels. A population inversion is obtained between the two 
lower ground state levels by submitting the Fe 3+ ions to a 
31 GHz signal. The maser oscillation takes place at 12.04 
GHz. As the sapphire resonator exhibits at low temperature 
extremely low losses, the Fe 3+ ion concentration required to 
get enough gain is very low. Contrary to the maser of the 60's, 
our WhiGMO (Whispering Gallery Maser Oscillator) incor- 
porates a microwave resonator presenting a Q-factor of typi- 
cally one billion at 4.2K made from a sapphire crystal contain- 
ing typically 0.1 ppm or less of active ions. The experimental 
validation of the 12 GHz WhiGMO and its potentiality to be 
competitive with the state-of-the-art microwave stable refer- 
ences opens the way to revisite the Maser concept and the 
associated material properties. Moreover the WhiGMO prin- 
ciple can be extented to other transitions of the Fe i+ ion in 
sapphire or to other doped crystals to design stable sources at 
higher frequencies. 



In this paper we refined the description of the Fe 3+ ion em- 
bedded into a (U-AI2O3) matrix to derive its 252 energy lev- 
els with an improved accuracy compared to previous works- . 
Moreover we calculated the Fe 3+ ion strength transition prob- 
ability including those for optical transitions. It appears than a 
few allowed optical transitions are favorable to obtain a maser 
operation at 3 1 GHz applying optical pumping at wavelength 
achievable with available low cost laser sources. Eventually, 
a first preliminary evaluation of the 3 1 GHz maser concept is 
presented. It is shown that the 31 GHz maser signal power 
achievable is larger than the classical 12 GHz version, i.e. a 
power gain of 100 is expected in the optimized conditions. As 
fundamental limit of frequency stability is fixed by the ther- 
mal noise, any increase in the maser signal output will benefit 
to the performance. 
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II. THE MOLECULAR HAMILTONIAN 

The electronic configuration of the Fe i+ cation is [Ar]3c/ 5 . 
The spectroscopic term of the ground state is 6 S (S= 5/2 and 
L= 0). In a doped sapphire crystal the Fe i+ ion replaces the 
cation Al i+ and undergoes the action of the crystal field of the 
nearby ions 2 ~ and AZ 3+ (figure[T). The cation is located in a 
distorted octahedral. Assuming the cation Fe i+ in symmetry 
Oh, the fundamental term will be 6 Ai s (<9/,)- 

Since the appearance of the fundamental articles of 
Racah 1 -, Eliot and Stevens^ and Judd^ 2 -, the theoretical inter- 
pretation of electronic spectra of transition metals has became 
increasingly easy. The appropriate molecular Hamiltonian is: 



+ Trees correction +Racah correction 



(1) 



In this section, we will discuss the electronic energy levels 
of a transition metal ion in crystalline solid. The theoretical 
foundation is already established and discussed in detail in 
the bibliograph y 12 " 16 . The Hamiltonian includes : the elec- 
trostatic Coulomb interactions (H ee ), the spin-orbit coupling 
(H so ), the ion-ligand interactions described in the framework 
of crystal field theory (H c ), the Zeeman interaction, and the 
Trees and Racah corrections. 



A. Free ion 

Assuming the radial function of the Fe 3+ d orbital as pro- 
posed by Zhaoii using the "double-zeta exponential" method, 
we find: 

R d {r) = 0.677 (^-)^ 2 r 2 exp(-5.6r) 



0.55237( 



6! 

3.446 7 , 



6! 



r exp(—\.132r) 



(2) 



From R c i(r), the Racah parameters Bo and Co, the spin-orbit 
coupling constant i^o, and the mean values {r k )o have been 
calculated for the free ion: 



So = 1130.22 cm- 1 , C =4111.45 cm" 1 



Co = 588. 946 cm 



-l 



(3) 



(r 2 ) = 1.89039 a.u. , (r 4 } = 11.46485 a.u. 

We take into account the Treesi£ and the Racahi^ correc- 
tions. The constant of Trees correction oco= 40 cm -1 and the 
Racah correction po= -131 cm~ l are fitted from the spectra of 
the free ioiUi 



B. Crystal field interaction 

According Wybourne 14 , the potential of the crystal field can 
be written as a linear combination of irreducible tensor oper- 



ators C„ 



(4) 



where B^q are the crystal field parameters. In the case of d- 
electrons, the allowed values for the integer k are 0, 2, 4 and 

-k<q< +k. 

The structure of the equation <j4j is determined by the sym- 
metry of the crystal lattice. For a trigonal field (C3,,), the sum- 
mation contains four terms: B20, B40, £43 and Z?4_3. In the 
generalized model of the crystal field, the trigonal crystal field 
parameters B^q were given by Zhao- The values of the re- 
duced matrix elements of tensor operators are published in the 
book of Nielson and Koster 2 - for 16 x 16 matrices. From the 
Wigner-Eckart theorem we developped our generalized model 
based on 252x252 matrices. 



C. Calculation method 

The general expression for the interaction between the cen- 
tral metal ion and ligand orbitals is complex, a reasonable ap- 
proximation for the electrostatic parameters B and C, the con- 
stant of spin-orbit, and the average values (r 4 ) was made that 



B=N 4 B , C=JV 4 C 
C=^V 2 Co , <r 4 )=jV 2 (r 4 )o 




a=N 4 a , p = JV 4 p 



(5) 



where N is the average reduction factor due to the covalency. 

In the AI2O3 crystal the position of Al i+ is defined respec- 
tively to the sites of the two nearest cations as schematized in 
the figure Q] 



a 




kq 



FIG. 1. Local structure distortion of octahedral Fe 3+ center in rx- 
A/2O3 system. 

The intersite distances and the corresponding angles with 
respect the symetry axis C3 have been determined in previous 
work s 21 ' 22 as: 



3 



R\=\.966A , 00=47.44° 

^=1.857i , 9^= 62.70° (6) 

To represent the crystal distortion due to the substitution 
of Al 3+ by Fe i+ , we assume as proposed by Zhao22 that the 
two distances is affected by an identical relative variation. We 
add an identical constraint to the angle relative variations, thus 
introducing two independent parameters / and g defined as 

f= R 1 = R 1 . = ^i = 5i 

J r>0 nO ' 8 „() fl o ( ' 

"1 ^2 1 2 

Eventually the three independent parameters N,f and g 
have been varied to approach the experimental observations. 
We found the optical spectrum of Fe 3+ in a — AI2O3 is well 
represented by taking : 

N = 0.863 , g = 0.996 and f = 1 .003 (8) 

The comparison between our calculations and experimental 
observations is shown in Table|I]compared to the values given 
by Zhao 9 , which the method serves as the basis of our own 
model. The introduction of a third independant adjustable pa- 
rameter leads to an improvement in the determination of al- 
most all the optical wavelengths in the Fe i+ ion spectrum. 

To strengthen our confidence in our physical description, 
we derived from our model the spin-Hamiltonian parameters 
describing the ground state 6 Ai(S) of the fe 3+ ion. With- 
out applied static magnetic field, there are three energy lev- 
els labelled in the following |/) with i = 1,2 or 3 (see figure 
121). From the expression of spin-Hamiltonian given in^i, the 
energy difference between the three hyperfine levels of the 
ground state are written as : 

£3 --Ei = ^[(18D + fl-F) 2 + 80a 2 ] 1/2 (9) 
E 2 -E\ = ^(fl--F)-D + i[(18D + a-£) 2 +80a 2 ] 1/2 
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TABLE I. Calculation and observed of the d-d transition for the crys- 
tal a-AhO?, : Fe 3+ (the values are in cm~ l ) . 
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FIG. 2. Structure of Fe i+ ground state in the OC-A/2O3 crystal without 
applied static magnetic field. 

The calculated values of the spin-Hamiltonien parameters: 
D, a, a-F and the energy gaps are given in Table HT1 and com- 
pared to the previous evaluation due to Zhao and experimental 
values. 
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TABLE II. Calculated spin-Hamiltonian parameters for OC-A/2O3 : 
Fe i+ at 300 K compared to experimental values (the unit is 10~ 4 
cm~' 
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III. TRANSITION PROBABILITY 

The optical spectrum of the iron doped sapphire shows 
fairly intense absorption bands 2 ^, 8 indicating that we can in- 
teract optically with the Fe 3+ ions. Our objective is to evalu- 
ate the possibility to get a population inversion in the ground 
state by applying an optical radiation to the crystal. Thus we 
need to calculate the effect of an optical radiation on the pop- 
ulation of each ground state level with i = 1,2 or 3 taking 
into account the allowed transitions with the different excited 
levels \j) listed in the table Q] Following Siegman- the transi- 
tion rate Wjj between one ground state level | ;} and one excited 
level \j) is given by: 

Wij = ^(ia) 2 g(v)a 2 (10) 

Where, y: the gyromagnetic factor. H: the module of the RF 
or optical magnetic field. g(v): the line shape function de- 
scribing the spectral broadening of the absorption line. 

G 2 is a dimensionless factor representing the strength of the 
transition probalility, which depends on the relative orienta- 
tion of the oscillating magnetic field with respect the crystal 
axis z- O 2 is computed from the matrix element of the in- 
teraction hamiltonian, which depends on the total magnetic 
momentum and on the applied oscillating magnetic field. In 
our case where no static magnetic field is present, two differ- 
ent orientations have to be considered: HizorH//z. We 
found that for each transition only one configuration gives no- 
ticeable value for a 2 . In the table |TTT] we give the relevant o? 
parameters for each transition. The relative orientation of H 
with respect z is indicated by the symbol (_L) or (//). 



IV. PROPOSAL FOR A 31 GHZ MASER 

Laser sources emitting some 100 mW at 555 nm with a 
linewidth less than 0. 1 nm are currently commercially avail- 
able (see for exi^L). From the table [TTT] we selected the tran- 
sition to the 4 E{G) excited level, at 18015 cm -1 (w 555 nm). 
The transition probabilities for the 6 Ai(S) — > 4 E(G) transi- 
tion with H_U are: a\= 1.40xl0~ 3 , o|= 1.80x 10~ 4 and rj 2 , 
« 0. Thus if we submit the crystal to a narrowband 555 nm 
laser source, the |1) level will be more effectively pumped, 
and thus large population differences between 1 1 ) and the two 
other ground state levels can be obtained. 



A. Rate equations 

Let'us consider a Fe i+ doped AI2O3 crystal submitted to 
a 555 nm narrowband source resonant with the transition 
6 A\(S) — > 4 E(G). The proposed system is schematized in 
figure [3] A cylindrical cryogenic sapphire resonator present- 
ing a high-Q resonance mode at 3 1 GHz is submitted to an 
optical radiation at 555 nm. The crystal axis z is colinear to 
the cylinder axis. The 31 GHz resonant mode is assumed to 
be quasi-transverse magnetic, i.e. the magnetic field lines in 
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TABLE III. Strength of transition probability a 2 for the Fe i+ d-d 
transitions in a — A/2O3. The integer (n) in brackets is for x 10~"). 



the equatorial plane of the cylinder are perpendicular to the z 
axis. A small magnetic loop placed near the dielectric cylinder 
is used to derive the 3 1 GHz maser signal output. 

The excited levels decay rapidly to the ground state with 
a lifetime To < 1 ms. All relaxation in the ground state are 
governed by the spin-lattice effect with a characteristic time 
Ti ss 10 ms at 4 K. It is also well established that the spin-to- 
spin relaxation is very fast with a characteristic time T2 < 20 
ns. In this case a simple rate equations model dealing only 
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FIG. 3. Energy levels and schematic of the cryogenic optically 
pumped 31 GHz WhiGMO . 



The total absorbed power is written as: 

Po = hVuWp An Uo NV e ff 
3 

= —hvuANuNVeff 



(13) 



N is the Fe i+ ions concentration and V e ff is the volume occu- 
pied by the active ions participating to the maser signal. V e /f 
is the volume of the 3 1 GHz mode in the sapphire resonator 
(see section HV Dl ). Here, we assume the optical radiation only 
lights the resonator part where the 3 1 GHz mode is concen- 
trated. Po is thus the minimal optical power required to get 
maser action in optimized conditions. 



with the level populations and neglecting coherences is suffi- 
cient to describe the interaction 30 . Excited ions can possibly 
fall on intermediate levels, nevertheless as To -C Xi they re- 
lax rapidly to ground state. This allows us to only consider a 
four levels system, i.e. the three ground state levels and the 
A E(G) excited one, which is labelled |4). The ground state 
ions submitted to the optical radiation will be excited with an 
absorption rate Wp, which is evaluated from equation[10] The 
optical pumping offers the possibility to get a population in- 
version between 1 3) and 1 1 } and thus to obtain a maser effect at 
3 1 GHz providing the resonator has been designed to present 
a high Q resonance at this frequency. The first 3 1 GHz emitted 
photons will stimulate others emissions with a rate Ws- If the 
optical pumping is effective enough to compensate for the res- 
onator loss, a 3 1 GHz self-sustained oscillation will occur. To 
theoretically determine the conditions to get a 3 1 GHz maser 
signal we followed the method described in 3 ! and validated in 
the case of 12 GHz WhiGMO. The stationary solutions of the 
rate equations are : 



A«i 3 = 



27AWi3 + 6(2AW 13 -3AWi 4 )xiW/> 
9(3 +4tiW s ) +4(3 +XiW s )xiWp 



A;; 
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9(3ANi4 + 2(2ANi4 - AM 3 )ti W s ) 
9(3+4tiW s )+4(3+XiW s )xiWp 



(Ha) 



(lib) 



ANjj and An,-,- are the population difference between levels 
|/) and \j) at the thermal equilibrium and in presence of the 
pump and maser signal respectively. The pump and maser 
signal frequencies are denoted V14 and Vi 3 respectively. 



B. Threshold pump power 

In the absence of the maser signal, the threshold pump 
power is achieved when the population inversion takes place. 
The absorption rate Wp corresponding to this situation is 
found by making A«i 3 = in the equation dl lab : 



Wp = 



AN 13 



2xi 3AN l4 -2AN 



(12) 



C. Maximum maser signal power 

A«i 3 and thus the 31 GHz signal power increase with the 
pump signal power system until saturation of the optical tran- 
sition occurs, i.e 114 ~ n\. Assuming we already overtook this 
situation, i.e. Wp — > °°, equations (11) simplify and the maxi- 
mum maser power is derived as: 



p 3W = hv 13 Max[W s Anu]NV e ff 
3 



(14) 



2x 



-/2Vi 3 (3ANi4 - 2ANn)NV e ff 
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Pj,\ max is the total power generated by the active ions inside 
the resonator. The coupling loop derive only a fraction of this 
power, which depends on the coupling coefficient p related to 
the loop area and to its position with respect to the sapphire 
resonator. For a critical coupling , i.e. (3 = 1, the maser output 
power is one half of P 3 i„, or . The resonator Q-factor is also 
impacted by the loading to the external circuit and for P = 1, 
the loaded Q-factor is one half of the unloaded one. 



D. Resonator design 

Our objective is to use the Fe 3+ ions present as contaminant 
in a high purity sapphire monocrystal. Typical effective ions 
concentration in such a type a crystal is well below 1 ppm. 
Thus the 3 1 GHz resonator has to present a very high quality 
factor, otherwise the 31 GHz stimulated emissions wont be 
sufficient to compensate for the resonator loss. This can be 
accomplished by using a whispering gallery mode, which 
Q-factor is only limited by the sapphire losses. Quasi-TM 
Whispering gallery modes are characterized by three integers 
m,n,l representing the electromagnetic field component 
variations along the azimutal, radial and axial directions 
respectively™ In the following we consider only the resonant 
modes with low radial and axial variations , i.e. those corre- 
sponding to n = I = 0. At low temperature, i.e. near the liquid 
helium temperature, the whispering gallery mode Q-factor 
can be as high as 1 billion providing the mode order is suffi- 
ciently high, i.e. m > 15 3 - 3 -. The WhiGMO principle has been 
demonstrated at 12 GHz, i.e. the frequency of the |1) — > |2) 



6 



transition. In this case, the population inversion results from 
a pumping at 3 1 GHz. The resonator has a diameter 2R = 50 
mm and a thichness /; = 30 mm. The quasi-TM Whispering 
Gallery mode WG£/n,o,o is used to support the maser signal^ 4 . 

For the whispering gallery mode the electromagnetic fields 
are concentrated near curved air-dielectric interface, between 
the resonator radius R and a caustic of radius Rc- The effective 
volume V e ff = %h (R 2 — Rc) occupied by the mode is a little 
mode dependant as Rc is given by: 



R C = 



/Vl3\ 2 



(15) 



2h 



where £,- « 9.4 is the sapphire permittivity is the transverse 
direction with respect z- c the light velocity in vacuum. 

To design a 3 1 GHz resonator we keep the same ratio di- 
ameter/thickness, which is near the optimal valued. Thus to 
design the cryogenic resonator only means choosing the az- 
imuthal number for the mode supporting the 31 GHz signal. 
m has to be higher than 15 to prevent radiation losses which 
limit the Q-factor and should not be to high otherwise the cou- 
pling with the output probe will be difficult to adjust. The ta- 
ble [IV] gives for different values of m, the resonator diameter 
2R, the threshold optical power Pq and the maximum maser 
power emitted by the ions. Here, the Fe 3+ concentration is 
assumed to be 0.02 ppm, i.e. the same value than the 12 GHz 
WhiGMO. 



m 


2R 


Po 






(mm) 


(mW) 


OuW) 


17 


19.22 


1.17 


1.02 


19 


21.13 


1.42 


1.24 


21 


23.05 


1.71 


1.49 


25 


26.88 


2.35 


2.05 


35 


36.46 


4.40 


3.85 


49 


50.00 


8.41 


7.36 



TABLE IV. Different resonator designs to support the 31 GHz signal. 
The pump threshold power Po and the maximum maser power Pji 
are given for few azimutal numbers as well as the resonator diameter 
assuming a constant ratio 2^ = | . 
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